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Minimum time generation of SU(2) transformations 
with asymmetric bounds on the control^ 
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We study how to generate in minimum time special unitary transformations for a two-level quantum 
system under the assumptions that: (i) the system is subject to a constant drift, (ii) its dynamics can 
be affected by three independent, bounded controls, (iii) the bounds on the controls are asymmetric, 
that is, the constraint on the control in the direction of the drift is independent of that on the 
controls in the orthogonal plane. Using techniques recently developed for the analysis of SU{2) 
transformations, we fully characterize the reachable sets of the system, and the optimal control 
strategies for any possible target transformation. 
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I. INTRODUCTION 

The implementation in minimum time of specific trans¬ 
formations is a key ingredient of many protocols requir¬ 
ing the manipulation of two-level quantum systems (as 
in quantum information theory [l|, in quantum optics, 
or in atomic and molecular physics). Usually, one is in¬ 
terested in mapping in minimum time an initial state to 
a final state under specific conditions (see e.g. and 
references therein). However, a general and convenient 
approach to this problem consists in considering as con¬ 
trol target the transformations themselves, rather then 
the states of the underlying physical system. In our con¬ 
text, the problem can be formulated as an optimal control 
problem on the Lie group S'C/(2) of special unitary trans¬ 
formations in 2 dimensions [5l-l7|. An arbitrary element 
of this group can be written as 


where a, (3 € C satisfy |ap -I- |/3p = 1. The dynamics of 
X is given by the Schrodinger operator equation 

X = [lUqJz + UxJx + '^yjy + UzJz)X, A(0) = I, ( 2 ) 

where Jk [k = x, y, z) are the skew-Hermitian generators 
of SU (2), that is, independent elements of the Lie algebra 
su( 2 ), (jjQ is an arbitrary real parameter characterizing a 
constant drift term in the dynamics, and Uk = Uk(t) are 
possibly time-dependent control actions constrained by 

2 I 2 ^ 2 2 ^ 2 /o\ 

Uy ^ 5 ^ 72 ) ( 3 ) 

where 71 ^ 0, 72 0. In other words, we assume that the 

strengths of the controls affecting the dynamics through 
Jz, or rather through generators depending on Jx and 
Jy, are independent. 
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Bounds on the controls depending on their squares, as 
in ®, naturally arise when the control actions are phys¬ 
ically realized through fields, with energy proportional 
to their square amplitude. The choice of bounds in (|3]) 
corresponds to different driving strengths along the z di¬ 
rection, or along directions in an orthogonal plane. This 
means that there is an anisotropy in the problem, with a 
privileged direction in space, determined by the specific 
apparatuses which are used to steer the system. Stronger 
constraints can be given by choosing independent bounds 
for the three controls, ^ uf, Uy ^ 7 I and ul < yf, 
and have been considered elsewhere (for instance, see 
0 and references therein). From a theoretical point of 
view, the analysis of the specific constraints considered 
in this work is relevant because it represents an interme¬ 
diate situation between the problem with independently 
constrained controls, which has not been solved in the 
general case, and the problem with the isotropic bound 
+ Uy + ul ^ 7 ^, which has been fully investigated 
[a, 0| • Physically, bounds on the controls depending on 
their squares are generically relevant in quantum infor¬ 
mation processing [ij, in atomic and molecular physics, 
and in Nuclear Magnetic Resonance (NMR) Q. 

The generators satisfy the standard SU(2) commuta¬ 
tion relations 

[Jj, Jk] = Ji, (4) 

where {j,k,l) is a cyclic permutation of (x,y,z). More 
explicitly, Jk = —^(Xk, where crfc are the Pauli matrices. 

The target of the control action is to steer the iden¬ 
tity I = A(0) to an arbitrary final operator Xf = X{tf) 
in minimum time tf, through a suitable optimal con¬ 
trol strategy Uk{t). To determine this strategy, we will 
use the necessary condition of optimality provided by the 
Pontryagin Maximum Principle (PMP) @,[13 , which we 
briefly review. We introduce an auxiliary variable, the 
so-called costate M £ su(2), represented by the coeffi¬ 
cients 

bk =-{M,X^JkX), k = x,y,z, (5) 

where bk = bk{t), and {A,B) = Tr(Ai?f). Then we 
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define the Pontryagin Hamiltonian as 

H{M,X,Vx,Vy,v^:) = ujobz + Vxbj; +Vyby + v^bz. ( 6 ) 

The PMP says that, if a control strategy Uk{t) satisfying 
the bounds ([ 2 ]) and the corresponding trajectory X(t) are 
optimal, then there exists a costate M 7 ^ 0 such that 

H(M, X, Ux, Uy,Uz) ^ H{M, X, Vx,Vy,Vz) (7) 

for all Vk satisfying ©. The PMP is only a necessary con¬ 
dition for optimality, useful for finding extremal control 
strategies and trajectories. The optimal strategy and tra¬ 
jectory are determined by comparing the extremal ones, 
analytically or numerically. 

For its relevance, the control of SU{2) operations has 
been extensively studied, with several constraints on the 
control protocols. In this work we follow the approach 
presented in [^, which allows an analytical investigation 
of the problem. In that paper, the cases with three or two 
controls were considered, with bounds u^ + Uy + ul ^ 7 ^ 
and + Uy < 7 ^ respectively. We refer to this work 
for several technical details which are omitted here for 
brevity. Note that the present framework reduces to the 
case with two controls in iia when 72 = 0. Moreover, 
some optimal solutions in the case with three controls 
in 0 are also optimal solutions in the present case, with 
values of 71 and 72 depending on the specific trajecto^. 
Therefore, our analysis complements that presented in [y, 
0 . 


II. DETERMINATION OF EXTREMAL 
TRAJECTORIES 


and Uz has values in the interval [— 72 , 72 ]- Finally, in 
case (iii) both ([5]) and (ITUl) must hold. 

Having the form of the extremal controls, we can 
integrate @ and determine the trajectories in SU{2). 
Case (i) is trivial, and we obtain X{t) = 
where t = |. We do not further consider this situa¬ 
tion which does not provide optimal solutions. Cases (ii) 
and (iii) can be jointly integrated, by remembering that 
—72 ^ U 2 ^ 72 or Uz = ±72 in the two cases, respec¬ 
tively. By using these controls in ([5]), it is possible to 
prove that 


bx = tLCOsiiot + (j)), by = yLS\n(u}t + (f), ( 11 ) 


where 4> is a constant, and w is a possibly time-dependent 
function given by 


1 f* 

= uj{t) = Wo + - / Uz{s) 
t JQ 


ds 


in case (ii), and the constant 


UJ = UJo + Uz - Jl- 




( 12 ) 


(13) 


in case (iii). Its range depends on the case under investi¬ 
gation: w < Wo -f 72 if > 0 , w > Wq — 72 if bz < 0 , and 
Wo — 72 ^ w(f) ^ Wo - 1-72 if bz = 0. We shall use w rather 
than bz and fj, to identify extremal trajectories. Notice 
that we have been able to integrate ([ 8 |) even in the case 
of time-dependent Uz (and then w) because of the simple 
form of this system. By substituting (ED into (nnD , we 
find 


Ux = cos {uit + (j)), Uj, = 7 i sin (wt-I-(/)). (14) 


The differential equations describing the costate dy¬ 
namics can be found by using ([4]) in ®, and they are 
given by 

bx = -{ujo + Uz)by + Uybz, 

by = {uJo + Uz)bx-Uxbz, ( 8 ) 

bz — Uxby Uybx- 

After defining fj, = ^b^ + b'^, we find that + b1 is a 

constant, which cannot vanish because M ^ 0. By max¬ 
imizing the Pontryagin Hamiltonian ([S]) with the con¬ 
straints ®, we determine the form of the extremal con¬ 
trols. There are three possible cases: (i) y = 0; (ii) 
bz = 0; (iii) neither of them. In case (i), it follows from 
(|8D that Ux = 0 and Uy = 0, and then bz is a non-zero 
constant. Therefore, it must be 

Uz = ^2siga{bz). (9) 

In case (ii), from maximization of H we determine the 
form of the extremal controls 

bx b,i 

Ux = 7l — , Uy = ll — , (10) 


By using these expressions in @, and considering the 
representation of X given in o, following the steps de¬ 
tailed in 0 (which can be simply readapted when w is a 
function of time), we find 

a = ^ cos or — i—sin ar^ 

(3 = —i—sinar, (15) 

a 

where we have rescaled time as r = |, and defined 

a = a{uj) = yjb'^ -kyj, (16) 

with 6 = 0 in case (ii), b = b{uj) = wq -f Uz — w in case 
(iii). 

Because of the form of the drift term and the bounds 
on the controls, the problem has a natural cylindrical 
symmetry. This is apparent from (ED, where the phase 
of /3 can be arbitrarily modified through the parameter 
(f). In other words, all operators X differing by an off- 
diagonal phase are completely equivalent in the frame¬ 
work adopted in this work, and they can be reached in 
the same optimal time 0. Consequently, we can fully 
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describe the extreme trajectories in SU (2) by consider¬ 
ing the evolution of a, or, more conveniently, its real and 
imaginary parts x and y respectively, which must satisfy 
x'^ + ^ 1- Therefore, we can represent the relevant 

trajectories in the unit disk of (or, equivalently, in 
C). When bz ^ 0 they are given by 

b , , 

a;± = cos WT cos ar-sinwrsinar (17) 

a 

b 

y± = — sin LOT cos ar - cos lot sm qt , 

a 

with u! as in (IH; when bz = 0 they are 

xo = coswr cosyiT, j/o = “ sinwr cosyir, (18) 
with Lo = Lo{t) as in (IT^ . 

In the analysis of extremal trajectories, it is often 
important to consider the so-called singular trajecto¬ 
ries, that is, extremal solutions such that the Pontryagin 
Hamiltonian is independent of the controls (for the rele¬ 
vance of extremal trajectories in concrete problems, see 
for instance [illil and references therein). By jointly 
considering m and dS]), we can conclude that, in the con¬ 
text considered in this work, these solutions do not exist, 
because they would require bx = by = bz = 0, which is 
inconsistent with the requirement M 0. Only regular 
trajectories (i.e., non singular) have to be taken into ac¬ 
count. Using a similar argument, we observe that it is 
impossible to concatenate the extremal trajectories de¬ 
scribed before. Again, this would require the vanishing 
of all the bj at the switching time, which is not admitted. 


III. EVOLUTION OF THE REACHABLE SETS 

Following [^, we define the optimal front-line as the 
set of terminal points for a candidate optimal trajectory 
at time r. As we have seen, depending on bz, there are 
three families of extremal trajectories. Correspondingly, 
there are three optimal front-lines in R^, 

^+{ t ) = {(x+,y+),-oo < w < Wo-f 72}, 

= {(a;-,2/-),wo - 72 < w < 00}, ( 19 ) 

^oi-r) = {(xo,?/o),wo - 72 < w ^ Wo-f 72}, 

or similar definitions in C, in terms of ao, a±. We remind 
that w is a constant for a given trajectory in J7|_ or J^_, a 
possibly time-dependent function for an extremal in J-q. 

The reachable set at time r is, by definition, the set of 
operators in SU{2) which can be reached in time smaller 
or equal than r. The evolution of the reachable set of 
the system is determined by the evolution of the opti¬ 
mal front lines (HU), in particular, by their intersections, 
where the trajectories could lose optimality. In the set¬ 
ting considered in Q , there is a unique optimal front line 
J^(t), and, if we work in C and forget for a while the 
bounds on w, J^±(t) can be expressed in terms of it as 

J-±(r) = e=F*'^=’'J'(T) ( 20 ) 


by means of suitable shifts in w. Therefore, assuming 
again that w can be any real number, we can write 

J-_(t) =e 2 *^^"J-+(r), ( 21 ) 


that is, at time r the two sets are mapped into each other 
by a rotation of angle 272 T in the unit disk. We observe 
that -F_|.(t) = when t = ^, with k &h. 

The individual analysis of J^+{t) and .F_(t) follows 
from that of J^(r). We summarize the main results. First 
of all, there is a one-to-one correspondence between val¬ 
ues of w and points on J7|_ and J--, that is, the associ¬ 
ated trajectories do not intersect in optimal conditions. 
Moreover, for each locus there is a critical trajectory [l^ 
spiraling around the center of the disk, modified with re¬ 
spect to that corresponding to J-'(r) according to ( 1 ^ . 
and parameterized by the critical frequencies 


7 i -f (wo ±72) 

<j-'c —-;-, 

Wo ±72 

and losing optimality at the critical times 


7 r|wo ±72! 

7i\/(wo ± 72 )^ +7i 


( 22 ) 


(23) 


These trajectories can be cut loci for the system, that 
is, special lines where optimal trajectories lose their op¬ 
timality. Other trajectories lose optimality on the border 
of the unit disk, which is, then, a cut locus for the sys¬ 
tem. At time r, the frequencies corresponding to these 
trajectories are given by 


Wc'( t) = (wo ± 72 ) ± -ll- (24) 

Note that, in (1^ and (1^ . quantities with sign -|- or — 
refer to or JA, respectively. The same applies for the 
first ± sign in (l24l) . but the second ± sign depends on 
the specific case. It is possible to prove that the possible 
scenarios are Wc > wq -f 72 > 0 or Wc < wq -f 72 <0 for 
J7|_, and Wc > wq — 72 > 0 or Wc < wq — 72 < 0 for T-. 
Therefore, considering the allowed range of values for w, 
we see that sometimes the critical trajectories are not 
extremal trajectories for the system. We can also refine 
the definition of the optimal front lines, neglecting con¬ 
tributions which are certainly sub-optimal. For instance, 
when Wq > 0 , the range of values of w for J^+(t) is given 
by Wc'(r) < LO < loq-\-^ 2 - For F-{t), it is Wc < w < Wc'(r) 
when 72 < Wo, and wq — 72 < w < Wc'(t) when 72 > wq. 
Similar expressions can be found when wq < 0. 

For small times (that is, in a neighborhood of t = 0) 
there is a one-to-one correspondence between values of 
w and points in Fq. Nonetheless, from (fT^ we see that 
there are different control strategies Uz = Uz{t) leading 
to the same w, that is, those having the same time aver¬ 
age. Therefore, in this case there are different trajecto¬ 
ries conver ging to the same point of Tq, and they are all 
equivalent |16l |. In other words, in the region spanned by 
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J-Q, there are distinct extremal trajectories (correspond¬ 
ing to different control strategies) leading to the same 
final state in the same time, and remaining extremals af¬ 
ter they intersect. This is not in contradiction with the 
results of where the optimal solution is unique and 
Uz is constant, because, in general, the optimal solutions 
for the case of asymmetric bounds parameterized by 71 
and 72 are not optimal solutions for the problem with 
symmetric bound given by 7^ = 7i + 72 ■ 

Since Xq -|- j/q = cos^ 71T, the optimal front-line 
is an arc of circle centered at the origin, with time- 
dependent radius COS71T, and angle 272T. If 2 ^ 2 T ^ 27 r, 
there are several values of w corresponding to the same 
point on the optimal front-line, and the corresponding 
extremal trajectories become equivalent, that is, they at¬ 
tain the same point at the same time time even if the 
time average of Uzit) is different. 

By considering the definition of and -^±(t), we 

see that these three loci are connected. Moreover, by 
using implicit differentiation, we find that 


dyo _ dy+ 
dxo dx+ 


dy- 

dx- 


cot LOT, 


( 25 ) 


therefore they are smoothly connected. It is possible to 
consider as optimal front line for this problem the union 
of these three loci. 

To complete the analysis, we must consider the inter¬ 
sections between T+(t), and J-q(t) at any time 

T. It turns out that J-q{t) never intersects T+{t) or 
T- (r) unless it is sub-optimal (and then these intersec¬ 
tions are irrelevant for the characterization of the evolu¬ 
tion of the reachable sets of the system). The intersection 
of -F+ (r) and T- (t) can be found numerically. The two 
endpoints of and , associated with co = coq + J 2 
and a; = wq —72 respectively, coincide in two cases: either 
when 'T = ^ (when the two endpoints of J-q overlap), or 
when r = 1^, when the radius of Tq vanishes. 


IV. TYPOLOGIES OF EVOLUTION OF THE 
REACHABLE SETS AND EXAMPLES 




FIG. 1: (Color online) Time evolution of the reachable sets 
in the unit disk, with 1^0 = 4, 71 = 1 (left plot) or 2 (right 
plot), 72 = 3. We have represented the optimal-front line 
at successive times t = 0.6, 1.0 and 1.4. The gray curve is 
JTi, the dashed and dotted lines represent the evolution of 
its endpoints, before and after they converge, respectively. 
The region spanned by Jh and enclosed in the dashed lines 
contains points which can be reached via several equivalent 
optimal protocols (and, possibly, with different time average 
of Uz (t) in the region enclosed in the dotted lines). The critical 
trajectory associated with T- is a cut locus for the system. 
The other cut loci, not shown in the plot, are the border of 
the unit disk and the set of intersections between J7|_ and J--. 


1/ 1/ 



FIG. 2: Time evolution of the reachable sets in the unit disk, 
with = 2, 71 = 1 (left plot) or 2 (right plot), 72 = 3, at 
successive times t = 0.6, 1.0 and 1.4. The cut loci for the 
system, not shown in the plot, are the border of the unit disk 
and the set of intersections between Jq. and J^_. 


We can sum up the previous results, and classify the 
systems in four classes, depending on the specific values 
of Wo, 7i and 72. They correspond to different forms 
of the optimal trajectories, producing different time- 
evolutions of the reachable sets. A given target operator 
Aj G SU{ 2 ) will require different control strategies (and 
associated minimum time tf) depending on the case at 
hand. For sake of simplicity we assume wq ^ 0 , but a 
completely analogous classification can be given also in 
the case wo < 0. 

First of all, we observe that, since wo -I- 72 ^ 0 the 
locus Jq. rotates counter-clockwise in the unit disk. The 
sense of rotation of J-- depends on the sign of wq — 72, 
therefore the evolution of the reachable set is radically 
different in the two cases 72 > Wq and 72 < Wq. Simi¬ 
larly, the long-time evolution of the reachable set depends 


on the relative magnitude of 71 and 72. Following the 
discussion of the previous section, if 271 ^ 72 there are 
extremal trajectories with ending points on Jq_ and J-- 
which can get arbitrarily close to the center of the unit 
disk (corresponding to the SWAP operator). Otherwise, 
when 271 < 72 there is a disk of radius costt^, centered 
about the origin, whose points can only be reached by 
trajectories associated with Jq. The relative magnitude 
between wo and 71 determines how many times the opti¬ 
mal front lines spiral around the origin before exhibiting 
the aforementioned features. For a graphical representa¬ 
tions of the evolution of the reachable set with several 
choices of the parameters see Fig.s [T] and [ 5 J In the first 
case we have 72 > wq, in the second case 72 < wq, with 
two different settings for 271 and 72 in both cases. 
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Note that, if 72 = 0 , the optimal front-line JFq collapses 
to a point and the analysis consistently reproduces that 
of @ in the case of two independent controls. 

To further illustrate the behavior of the reachable sets, 
we provide some examples of optimal synthesis of some 
standard unitary gates, and compare with the corre¬ 
sponding results in the case of symmetric bounds. 

In the case of diagonal target operators, Xf = 
with A G [ 0 , 27 r), since these are represented by points 
on the unit circle, which are reached by the trajecto¬ 
ries forming or , the optimal control strategies are 
given by controls Ux and Uy as in ( 1141 ) with uj = ujc' as in 
(l 24 l) . and Uz = ±72, for J-± respectively. From (l 20 l) and 
the results in Q, the optimal time is given by 




2 min 

li2=±72 


/(tt - A)(wo + Uz ) + 
V (a;o + Uz)^+7i 


( 26 ) 


where 7r'^(uJo + UzY + A( 27 r — A)7j. In particular, 

according to the former discussion on the evolution of the 
reachable sets, the minimum is obtained with Uz = 72 
when 72 < Wo; when 72 > wq the situation is more 
complicated. The optimal time ( 1 ^ and the correspond¬ 
ing strategies can be compared with analogous time and 
strategies in the case of a symmetric bound on the con¬ 
trols 0, 


7 


4-n--2A 

2A 

7—Wo ’ 


if Wo > 2 L ^7 
if Wo < ^7 


( 27 ) 


obtained with Ux = Uy = 0 and = 7 or —7. If we 
require 7^ = 7 i + 7! (that is, the total control strength 
is the same), we have that the time (1221) is smaller than 
(USD, since the scenario with asymmetric control bounds 
is compatible with a symmetric bound. As a check of con¬ 
sistency, this result can be proven by using the Lagrange 
multipliers method on (1^ . leading to the constrained 
minimum (EH) obtained when 71 = 0 and 72 = 7. 

As a second example, we choose as target operator the 
SWAP operator Xf = ia-y, which represents the NOT 
operation in quantum information. In this case, the be¬ 
havior of the optimal trajectories is described by the opti¬ 
mal front line J-,, which we have analyzed in the previous 
section. The optimal strategies are given by Ux and Uy 
as in () 14 l) with w taking an arbitrary value in the interval 
[wo — 72, Wo + 72], possibly time-dependent. The control 
Uz can take any form, in particular Uz = 0 can be cho¬ 
sen, leading to w = Wq. For any choice of the parameters, 
the SWAP operator is attained in optimal time tf = 
which is independent on 72, consistently with the fact 
that Uz is completely irrelevant for the optimal synthesis 
of this operator. In the case of a symmetric bound on the 
controls, the optimal control strategy has a similar struc¬ 
ture Q, with Ux = 7 cos (wot-b <p), Uy = 7sin (wpt-b <p) 
and = 0 (tp is a phase), and the optimal time is tf = 
Again, this optimal time is smaller than that with asym¬ 
metric bound on the controls, under the assumption that 
7^ = 7i + 72- It is clear that the case with symmetric 
bound is reproduced when 72 = 0 and 71 = 7. 


V. CONCLUSIONS 


We have fully characterized the evolution of the reach¬ 
able sets of system @ with controls subject to the asym¬ 
metric constraints We have derived the optimal 

trajectories and the optimal control strategies for the 
system, and provided its classification in terms of the 
dynamical parameters, which are completely arbitrary: 
driftless dynamics or unbounded control actions are spe¬ 
cial cases of this treatment. By using this analysis, it is 
possible to compute, at least numerically, the minimum 
time for generating an arbitrary SU{ 2 ) transformation, 
and the required strategy. For sake of clarity, we have an¬ 
alyzed the cases of diagonal operations, and the SWAP 
transformation of a two-level system. These examples 
clearly illustrate the role of the constraints on the con¬ 
trols in the study of time-optimal synthesis of quantum 
operations. 

The main tool for studying the evolution of the reach¬ 
able sets is the optimal front-line, which has already 
proven useful for the investigation of the minimum-time 
synthesis of SU{ 2 ) operations. Its application to similar 
problems on different Lie groups is of great potential in¬ 
terest. It does not only provide a way to clearly visualize 
the behavior of the reachable sets, but also a simple ap¬ 
proach to prove rigorous results, whose derivation when 
following the separate trajectories could be cumbersome 
in some regions of the space of dynamical parameters. 

For instance, the results in SU{ 2 ) can be used to per¬ 
form a similar analysis in SO{ 3 ) (because of the stan¬ 
dard homomorphism connecting these groups) and are 
therefore related to the problem of attitude control of a 
rigid body. In this context, asymmetric bounds on the 
controls, as those considered in this work, are especially 
relevant for the treatment of rigid bodies with rotational 
symmetry about one axis. 

Another problem where the specific investigations pre¬ 
sented in this paper could be of relevance is the synthesis 
of SU{ 2 ) operations with individual bounds on Ux, Uy 
and Uz- Also in this case the problem has more degrees 
of freedom, since the cylindrical symmetry of this work 
and [ 1 , 0 ] is broken. However, the analysis of suitable 
optimal front lines could provide new insights for the in¬ 
vestigation of regions, in space of dynamical parameters, 
which have not been studied so far. 

Generalizations of this technique to problems charac¬ 
terized by an higher number of degrees of freedom seems 
a promising research line. Despite in these cases it seems 
difficult to obtain the particularly simple representation 
of the evolution of the reachable sets arising in SU{ 2 ), 
we believe that an approach based on the study of the 
envelopes of the front lines is more promising than a di¬ 
rect analysis of the trajectories. A prospective direction 
for this line of research is the study of the optimal im¬ 
plementation of two-qubit gates, or, more gen erally, the 
simultaneous control of two spins (see [1^ [Ij| for some 
recent applications of the PMP principle in this context). 
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